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' In this article, we investigate the under-structures of the tt and B mesons in the framework 

of the Bethe-Salpeter equation with the confining effective potential (infrared modified flat bottom 
potential). In bare quark-gluon vertex approximation, we obtain the algebraic expressions for the 
solutions of the coupled rain-bow Schwinger-Dyson equation and ladder Bethe-Salpeter equation. 
Firstly, we neglect the rain-bow Schwinger-Dyson equation, take the bare quark propagator and 
solve the Bethe-Salpeter equation numerically alone. Although the bare quark propagator can not 
embody dynamical chiral symmetry breaking and has a mass pole in the time-like region, it can 
■ give reasonable results for the values of decay constants and fs compared with the values of 

^ ' experimental data and other theoretical calculations, such as lattice simulations and QCD sum 

, rules. Secondly, we explore those mesons within the framework of the coupled rain-bow Schwinger- 

' Dyson equation and ladder Bethe-Salpeter equation. The Schwinger-Dyson functions for the u 

and d quarks are greatly renormalized at small momentum region and the curves are steep at about 
— IGeV^ which indicates an explicitly dynamical symmetry breaking. The Euclidean time fourier 
transformed quark propagator has no mass poles in the time-like region which naturally implements 
confinement. As for the b quark, the current mass is very large, the renormalization is more tender, 
however, mass pole in the time- like region is also absent. The Bethe-Salpeter wavefunctions for both 
' the TT and B mesons have the same type (Gaussian type) momentum dependence as the corresponding 

C^, wavefunctions with the bare quark propagator, however, the quantitative values are changed and 

^ the values for the decay constants and fs axe changed correspondingly. 
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I. INTRODUCTION 



Quantum chromodynamics (QCD) is the appropriate theory for describing the strong interaction at high energy 
region, however, the strong gauge coupling at low energy destroys the perturbative expansion method. The physicists 
propose many nonperturbative approaches to deal with the long distance properties of QCD, such as Chiral pertur- 
bation theory [1], heavy quark effective theory [2], QCD sum rule [3], lattice QCD [4], perturbative QCD [5], coupled 



*E-mail,wangzgyiti@yahoo. com.cn 
t Mailing address 



1 



Schwinger-Dyson equation (SDE) and Bethe-Salpeter equation (BSE) method [6], etc. All of those approaches have 
both outstanding advantages and obvious shortcomings. For example, lattice simulations are rigorous in view of QCD, 
they suffer from lattice artifacts and uncertainties, such as Gribov copies, boundary conditions and so on, further- 
more current technique can not give reliable result below 1 GeV, where the most interesting and novel behavior is 
expected to lie. The coupled SDE and BSE have given a lot of successful descriptions of the long distance properties 
of strong interactions and the QCD vacuum, for a recent review one can see Ref. [7]. The SDE can provide a natural 
way to embody the dynamical symmetry breaking and confinement which are two crucial features of QCD, although 
they correspond to two very different energy scales [8,9]. On the other hand, the BSE is a conventional approach 
in dealing with the two body relativistic bound state problems [10]. From the solutions of the BSE, we can obtain 
useful information about the under-structure of the hadrons and thus obtain powerful tests for the quark theory of 
the mesons. However, the main drawback can be traced back to the fact that when we solving the SDE and BSE, 
model dependent kernels for the gluon two point Green's function have to be used, furthermore, the coupled SDE and 
BSE are a divergent series of equations, we have to make truncations in one or the other ways. Numerical calculations 
indicate that the coupled rainbow SDE and ladder BSE with phenomenological potential models can give satisfactory 
results. The usually used effective potentials are confining Dirac 6 function potential, Gaussian distribution potential 
and flat bottom potential (FBP) [11-13]. The FBP is a sum of Yukawa potentials, which not only satisfies gauge 
invariance, chiral invariance and fully relativistic covariance, but also suppresses the singular point which the Yukawa 
potential has. It works well in understanding the dynamical chiral symmetry braking, confinement and the QCD 
vacuum as well as the meson structures, such as electromagnetic form factor, radius, decay constant [14,15]. 

The decay constant of the B meson /b plays an important role in modern physics with the assumption of current- 
meson duality. The precise knowledge of the value of the fs will provide great improvement in our understanding 
of various processes convolving the B meson decays. At present, it is a great challenge to extract the value of the 
B meson decay constant fs from experimental data. So it is interesting to combine the those successful potentials 
within the framework of coupled SDE and BSE to calculate the decay constants of both the tt and B mesons. In this 
article, we use an infrared modified flat-bottom potential (IMFBP) which takes the advantages of both the Gaussian 
distribution potential and the FBP to calculate both the tt and B mesons decay constants. Certainly, our potential 
model can be used to investigate the properties of other pseudoscalar mesons, such as K,D,Ds - ■ ■■ For example, 
we can obtain the decay constants = 127MeV, fx = 156MeV, fn = 238MeV and /s = 192MeV with the 
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same parameters, while a detailed studies of those mesons K,D,Ds - ■ ■ may be our next work, they are not our main 
concern in this article. 

The article is arranged as follows: we introduce the infrared modified flat bottom potential in section II; in section 
III and IV, we solve the Schwinger-Dyson equation and the Bethe-Salpeter equation and obtain the decay constants 
for both the tt and B mesons; section V is reserved for conclusion and discussion. 



The infrared structure of the gluon propagator has important implication for the quark confinement. One might 
expect that the behavior of the quark interaction in the region of small space- like determines the long range proper- 
ties of the qq potential and hence implements confinement, however, the present techniques in QCD manipulation can 
not give satisfactory small r behavior for the gluon propagator, on the other hand, the phenomenological confining 
potential models give a lot of successes in dealing with the low energy hadron physics, such as dynamical chiral sym- 
metry braking, pseudoscalar mesons electromagnetic form factors, mass formulations, w — w scattering parameters, 
etc [7,11,16]. In this article, we use a gaussian distribution function to represent the infrared behavior of the gluon 
propagator. 



which determines the quark-quark interaction through a strength parameter zu and a ranger parameter A i. This 
form is inspired by the 6 function potential (in other words the infrared dominated potential) used in Ref. [11], which 
it approaches in the limit A — > 0. For the intermediate momentum, we take the FBP as the best approximation and 
neglect the large momentum contributions from the perturbative QCD calculations as the coupling constant at high 
energy is very small. The FBP is a sum of Yukawa potentials which is an analogous to the exchange of a series of 
particles and ghosts with different masses (Euclidean Form), 



where N stands for the minimum value of the masses, p is their mass difference, and aj is their relative coupling 
constant. 



Here we correct a writing error in the first version. 



II. INFRARED MODIFIED FLAT BOTTOM POTENTIAL 



47rG(fc2)=37r2^e-^ 



(1) 



n 




(2) 
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The definition of momentum regions between infrared and intermediate momentum is about Aqcd = 200MeV, 
which is naturally set up by the minimum value of the masses N = IAqcd, where the gaussian function e~"S" 
decays to about 0.3 of its original values. Certainly, there are some overlaps between those regions, in this way, 
we can guarantee the continuity of the momentum. The asymptotic freedom tell us that at high energy the gauge 
coupling is very small and can be neglected safely, on the other hand, our phenomenological potential at energy about 
N + jp,j > 3 is already extend to the perturbative region and catches some perturbative physical effects. Thus, our 
phenomenological infrared modified FBP is supposed to embody a great deal of physical information about all the 
momentum regions. 

Due to the particular condition we take for the FBP, there is no divergence in solving the SDE. In its three 
dimensional form, the FBP takes the following form: 

" e-{N+jp)r 

V{r) = -^a, . (3) 

In order to suppress the singular point at r = 0, we take the following conditions: 

V{0) = constant, 

dr dr"^ dr'^ 

So we can determine Uj by solve the following equations, inferred from the flat bottom condition Eq.(4), 

n 

aj = 0, 

3=0 

n 

3=0 

n 

3=0 



Ya^{N+jpr = 0. (5) 
3=0 



As in previous literature [13-15], n is set to be 9. 

III. SCHWINGER-DYSON EQUATION 



The Schwinger-Dyson equation, in eff'ect the functional Euler-Lagrange equation of the quantum field theory, 
provides a natural framework for investigating the nonperturbative properties of the quark and gluon Green's functions. 

4 



By studying the evolution behavior and analj^ic structure of the dressed quark propagator, one can obtain valuable 
information about the dynamical symmetry breaking phenomenon and confinement. The SDE for the quark takes 
the following form: 

1 Rtt-; r H'^k 

S-\p) = h-p + m+^ j ^r^5(fc)7.G^.(fc - P), (6) 

where 

S-'ip) = iAip^h ■ P + B{p^) = A{p%i^ ■ p + m(/)], (7) 
G,.{k) = {5,. - ^)G{k% (8) 

and m stands for an explicit quark mass-breaking term. In the rainbow approximation, we take = 'Jij,. With the 
explicit small mass term for the u and d quarks, we can preclude the zero solution for the B{p) and in fact there indeed 
exists a small bare current quark mass. In this article, we take Landau gauge. This dressing comprises the notation 
of constituent quark by providing a mass m{p'^) = B{p'^)/A{p'^), which is corresponding to the dynamical symmetry 
breaking. Because the form of the gluon propagator G{p) in the infrared region can not be exactly inferred from the 
SU{3) color gauge theory, one often uses model dependent forms as input parameters in the previous studies of the 
rainbow SDE [6,7,13-16], in this article we use the infrared modified FBP to substitute for the gluon propagator. 

In this article, we assume that a Wick rotation to Euclidean variables is allowed, and perform a rotation analyti- 
cally continuing p and k into the Euclidean region where them can be denoted by p and k, respectively. Alternatively, 
one can derive the SDE from the Euclidean path-integral formulation of the theory, thus avoiding possible difiiculties 
in performing the Wick rotation [17] . As far as only numerical results are concerned, the two procedures are equal. 
In fact, the analytic structure of quark propagator has interesting information about confinement, we will go to this 
topic again in the third subsection of section 4. 

The Euclidean rain-bow SDE can be projected into two coupled integral equations for A{p^) and B{p'^), the 
explicit expressions for those equations can be found in Ref. [14,15]. For simplicity, we ignore the bar on p and k in 
the following notations. 

IV. BETHE-SALPETER EQUATION 

The BSE is a conventional approach in dealing with the two body relativistic bound state problems [10]. The 
quark theory of the mesons indicate that the mesons are quark and anti-quark bound states. The precise knowledge 
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about the quark structures of the mesons will result in better understanding of their properties. In the following, we 
write down the BSE for the pseudoscalar mesons, 

S+\q + ^P)x{Q,P)Sz\q-{l-m = ^ j j^V ^x{k,P)V.G ^.{q - k), (9) 

where S{q) is the quark propagator function, G^vik) is the gluon propagator, is the four-momentum of the center 
of mass of the pseudoscalar mesons, is the relative four-momentum between the quark and anti-quark in the 
pseudoscalar mesons, is the full vertex of quark-gluon, ^ is the center of mass parameter which can be chosen to 
between and 1, and x{q^P) is the Bethe-Salpeter wavefunction (BSW) of the bound state. In the limit = 'Jn, 
we obtain the ladder BSE. 

After we perform the Wick rotation analytically and continue q and k into the Euclidean region, the Euclidean 
pseudoscalar BSW xiq^ P) can be expanded in lorentz-invariant functions: 

Xiq, P) = 75 [iFi{q, q-P)+j- PF^iq, q-P)+^- qFaiq, q ■ P) + ■ q,^ ■ P]Fi{q, q-P)]. (10) 

The BSW Fi can be expressed in terms of the 50(4) eigenf unctions, the Tchebychev polynomials T-i (cos^), 

oo 

Fi{q,q- P) = ^Fr{q,P)q^P^T^ {cos 0), (11) 



where n = even if i = 1,2,4; n = odd if n = 3, Tn {cos 9) = cos(ncos^) and 6 is the included angle between q and 
P. In solving the coupled BSEs for F", it is impossible to solve an infinite series of coupled equations, we have to 
make truncations in one or the other ways in practical manipulations, numerical calculations indicate that taking only 
n = 0, 1 terms can give satisfactory results, 

X{q, P) = 75 [^F^{q, P)+l- PF^{q, P) + 1 ■ qq ■ PFi{q, P) + i[7 • g, 7 • P^iq, P)] • (12) 

For a thorough investigation of the solutions of the above BSWs, we must take full quark propagator and full quark- 
gluon vertex, again we are led to solve a divergent series of coupled SDEs and BSEs, truncations in one or the other 
ways for the quark propagator and quark-gluon vertex are also necessary. In this article, we take the bare vertex for 
both the SDE and BSE. 

In solving the BSEs, it is important to translate the wavefunctions into the same dimension, 

pn _^ f^2npn^ pn _^ ^^n+l pn ^ pn _^ f^2n+l pn ^ pn _^ ^2n+2 pn ^ ^ _^ p _^ pj^ ^ (^3^ 
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where A is some quantity of the dimension of mass. 

Here we take a short digression to discussing the spectrum of the BSEs. In ideal conditions, a precise solution 
of the BSE for the bound states of definite quantum numbers will reproduce the full spectrum with the fundamental 
parameters of QCD, such as SU{3) gauge invariance, quark masses, etc. For example, the solutions of the BSE 
for mesons will result in a full pseudoscalar spectrum for both the fundamental states and excited states such 
as 7r°, 7r(1300) • • •. However, the present conditions are far from the case, the truncated BSEs always result in a 
spectrum with more bound states (artifact) [18]. Moreover, the spectrum is not the major subject which the present 
article concern. So in the article, we take the masses of the pseudoscalar mesons as input parameters and make an 
investigation of the ir and B mesons BSWs for both ladder approximation and bare quark propagator approximation. 

The ladder BSE can be projected into four coupled integral equations in the following form, 

H{l,l)F^{q,P) + H{l,2)F^{q,P) + H{l,3)F^{q,P) + H{l,4)F°{q,P)= / Pdk sm^eK{l,l), 

Jo Jo 

roo rir 

H(2,l)F°{q,P) + H{2,2)F°{q,P) + H{2,3)Fi{q,P)+H{2,4)F°{q,P)= / k^'dk sin'' 6{K(2,2) + K{2,3)), 

Jo Jo 

H{3,l)F^{q,P) + H{3,2)F^{q,P) + H{3,3)F^{q,P) + H{3,4)F°{q,P)= / k'^dk sin' e{K{3,2) + K{3,3)), 

Jo Jo 

roo rir 

H{A,l)F°{q,P) + H{A,2)F^{q,P) + H{A,3)F^{q,P)+H{4,4)F°{q,P)= k^dk sinHK{A,4), (14) 

Jo Jo 

the expressions of the H{i,j) and K{i,j) are cumbersome and neglected here, the interested readers can get the 
word- version from the author. 

Here we give some explanations about the expressions of H{i,j) . The H{i,jys are functions of the quark's 
Schwinger-Dyson functions (SDF) A{q' + ^'P' + ^q ■ P), B{q' + ^'P' + ^q ■ P), A{q' + (1 - ^^^2 _ _ . 
and B{q' + (1 — ^)'P' — (1 — • P). The relative four-momentum q is a quantity in Euclidean space-time while the 
center of mass four-momentum P is a quantity in Minkowski space-time. The present theoretical techniques can not 
solve the SDE in Minkowski space-time, we have to expand A and B in terms of Taylor series oiq- P, 

A{q' + C'-P' +iq-P)= A{q' + ^'-P') + A{q' + C'P')'^9 • P + • • • , 



B{q' + C2p2 +^q.p)= B{q' + ^P^) + B{q' + fP'Y^q • P + • • • . (15) 

The other problem is that we can not solve the SDE in the time-like region as the two point gluon Green's function 
can not be exactly inferred from the SU{3) color gauge theory even in the low energy space-like region. In practical 
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manipulations, we can extract the values of A and B from the space-like region smoothly to the time-like region with 
the polynomial functions. To avoid possible violation with confinement in sense of the appearance of pole masses 
= —m{q^), we must be care in the choice of polynomial functions [11]. 
Finally, we write down the normalization condition for the BSW, 



/ 



where x = 74%"^ 74- We can substitute the expressions of the BSWs and SDFs into the above equation and obtain the 
precise result, however, the expressions are cumbersome and neglected here. 

A. Decay constants of pseudoscalar mesons 

The decay constants of the pseudoscalar mesons are defined by the following current-meson duality, 

lUPi. = (0|g7M75<?k(P)), 



= ^JN, / Tr 



(2^ 



(17) 



here we use tt to represent the pseudoscalar mesons. 

B. Bethe-Salpeter equation with bare quark-gluon vertex and bare queirk propagator 

In this subsection, we investigate the BSWs of the tt and B mesons with the quark-gluon vertex and quark 
propagator are both taken to be bare, 

r^ = 7M> S-\p) = i^-p + M, (18) 

where the effective mass M is taken to be the constituent quark mass for the u, d and b quarks. In this article, we 
take the effective mass M as an input parameter. Such a two point quark Green's function can not embody dynamical 
chiral symmetry breaking and has a mass pole in the time-like region. However, as a first step, we can study the 
BSWs for those pseudoscalar mesons with the crude approximation. The algebraic expressions for the BSWs can be 
obtained easily with a simple substitution of A{p) 1 and B{p) ^ M in Eq.(14). After solving the BSEs numerically 
by iterations, we plot the BSWs Ff, , F^, as functions of the relative four- momentum q for the tt meson and B 
meson, respectively. In this article, we take the tt meson BSWs explicitly shown in Fig.l as an example and neglect 
others for simplicity. As the values of the wavefunctions F^ , F^ are tiny, we plot them perspicuously in another figure 



. From those figures, we can see that the first two wavefunctions are dominating, -F°, -F2' ^ -^3 ) -^4 > and all of the four 
wavefunctions are gaussian-type and center around very small momentum i.e. near zero momentum which indicates 
that the bound states must exist at the small momentum region or in other words confinement occurs at the infrared 
region. Based on the numerical values of the BSWs of the tt and B mesons, we can obtain the corresponding decay 
constants. 

U = IMMeV; fs = 164MeV, (19) 

which are compatible with the experimental, lattice and QCD sum rule results, = 130MeV{Exp) and fs ~ 
150- 210 MeV{Latt, sumrule) [20-22]. In calculation, the input parameters are N = l.OA, V(0) = -10. OA, p = 5. OA, 
M„ = Md = 530MeV, Mb = 5200Mey, A = 200Mey, w = 1.3GeV and A = O.OSGeV^. 

C. Coupled rain-bow Schwinger-Dyson equation and ladder Bethe-Salpeter equation 

In this subsection, we explore the coupled equations of the rain-bow SDE and ladder BSE with the bare quark- 
gluon vertex for both the tt and B mesons. The algebraic expressions for those solutions are obtained already in 
section 3 and forepart of section 4, here we will not repeat the tedious routine. In solving those equations numerically, 
the simultaneous iterations converge quickly to an unique value independent of the choice of initial wavefunctions. 
The final results for the SDFs and BSWs are plotted as functions of the square momentum q^. 

The quark-gluon vertex can be dressed through the solutions of the Ward-Takahashi identity or Salvnov- Taylor 
identity and taken to be the Ball-Chiu vertex and Curtis- Pennington vertex [23,24]. Although it is possible to solve 
the SDE with the dressed vertex, our analytical results indicate that the expressions for the BSEs with the dressed 
vertex are cumbersome and not suitable for numerical iterations 

In order to demonstrate the confinement of quark, we have to study the SDF of the quark and prove that there 
no poles on the real timelike axial. So it is necessary to perform an analytic continuation of the dressed quark 
propagator from Euclidean space into Minkowski space p4 ipo. However, we have no knowledge of the singularity 
structure of quark propagator in the whole complex plane. One can take an alternative safe procedure and stay 
completely in Euclidean space avoiding analytic continuations of the dressed propagators [25]. It is sufficient to take 



This observation is based on the authors's work in USTC. 
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the Fourier transform with respect to the EucUdean time T for the scalar part 5, 




ig^T q 



2n 



s ) 



(20) 



If S(p) had a pole at = — m^, the Fourier transformed S*{T) would fall off as e~™^ for large T or log 5* = —mT. 

In our calculation, for large T, the values of S* is negative, except occasionally a very small fraction positive 
values. We can express S* as |5*|e'"'^, n is an odd integer, log 5* = log IS'*! + inn. If we neglect the imaginary 
part, we find that when the Euclidean time T is large, there indeed exists a crudely approximated (almost flat) linear 
function with about zero slope for all the u, d (The curve for the d quark has the same behavior as the u quark in the 
limit of Isospin symmetry is exact ) and b quarks with respect to T, which is shown in Fig. 2. Here the word 'crudely' 
should be understand in the linearly fitted sense, to be exact, there is no linear function. However, such fitted linear 
functions are hard to acquire physical explanation and the negative values for S* indicate an explicit violation of the 
axiom of reflection positivity [26], in other words, the quarks are not physical observable i.e. confinement. 

From Fig. 3, we can see that for the u and d quarks, the SDFs arc greatly renormalized at small momentum region 
and the curves are steep at about = IGeV"^ which indicates an explicit dynamical chiral symmetry breaking, while 
at large g^, they take asymptotic behavior. As for the h quark, shown in Fig.4, the current mass is very large, the 
renormalization is more tender, however, mass pole in the time- like region is also absent, which can be seen from 
Fig.2. At zero momentum, rn„(0) = 688Mey, rnd(O) = 688MeT/ and rribiQ) = AdQOMeV, which are compatible with 
the constituent quark masses. In fact, the connection of m(p) to constituent masses is somewhat less direct for the 
light quarks and is precise only for the heavy quarks. For heavy quarks, niconstituentip) = m{p = 2mconsUtuent{p)) , 
for light quarks , it only makes a crude estimation [19]. From the plotted BSWs (neglected here for simplicity), we 
can see that the BSWs for both the tt and B mesons have the same type momentum dependence as the corresponding 
wavefunctions with the bare quark propagators, however, the quantitative values are changed. The gaussian type BS 
wavefunctions which center around small momentum indicate that the bound states exist only in the infrared region, 
in other words confinement. Finally we obtain the values for the decay constants of the tt and B mesons, 



which are compatible with the experimental, lattice and QCD sum rule results, f„ = 130MeV{Exp) and fs ~ 



127MeV; 



/b = 192MeV, 



(21) 
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150-210MeV{LaU, sumrule) [20-22]. In calculation, the input parameters are N = l.OA, V{0) = -ll.OA, p = 5.0A, 
m„ = ma = 6MeV, rrib = AlOOMeV, A = 200MeV, w = 1.6GeV and A = OMGeV^. 

From the variations of the values for the decay constants of both the n and B mesons, we can estimate that the 
full vertex approximation will not change those values greatly. 

V. CONCLUSION AND DISCUSSION 

In this article, we investigate the under-structures of the tt and B mesons in the framework of the Bethe-Salpeter 
equation with the confining eff'ective potential (infrared modified flat bottom potential). In bare quark-gluon vertex 
approximation, we obtain the algebraic expressions for the solutions of the coupled rain-bow SDE and ladder BSE 
for those mesons. At the first step, we neglect the rain-bow SDE, take the bare quark propagator and solve the BSE 
numerically alone. Although the bare quark propagator can not embody dynamical chiral symmetry breaking and 
has a mass pole in the time- like region, it can give reasonable results for the values of the decay constants f^^, Jb 
compared with the values of the experimental data and other theoretical calculations, such as lattice simulations and 
QCD sum rules. In calculation, we obtain the BSWs for the tt and B mesons, which center in the small momentum 
region, are compatible with confinement. Secondly, we explore the coupled equations of the rain-bow SDE and 
ladder BSE with the bare quark-gluon vertex for those mesons. The quark-gluon vertex can be dressed through 
the solutions of the Ward-Takahashi identity or Salvnov- Taylor identity and taken to be the Ball-Chiu vertex and 
Curtis-Pennington vertex, however, a consistently numerical manipulation is unpractical. After we solve the coupled 
rain-bow SDE and ladder BSE numerically, we obtain both the SDFs and BSWs for both the tt and B mesons. The 
SDFs for the u and d quarks are greatly renormalized at small momentum region and the curves are steep at about 
= IGeV^ which indicates an explicitly dynamical chiral symmetry breaking. After we take Euclidean time fourier 
transformation about the quark propagator, we can find that there is no mass pole in the time-like region and obtain 
satisfactory result about confinement. As for the b quark, the current mass is very large, the renormalization is more 
tender, however, mass pole in the time-like region is also absent. The BSWs for both the tt and B mesons have the 
same type momentum dependence as the corresponding wavefunctions with the bare quark propagators, however, the 
quantitative values are changed and the corresponding values for the decay constants /tt and /b are changed, but not 
greatly. We can estimate that the full vertex approximation will not change those values greatly. Once the SDFs and 
BSWs for both the tt and B mesons are known, we can use them to investigate a lot of important quantities in the B 
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meson decays, such as B — tt, B — K,B — p former factors, Isugar-Wise functions, strong coupling constants, etc. 
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Fig.1 BSWs of n meson 
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Fig.3 SDFs of u and d quark 
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